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Abstract 

Assuming the QCD multipole expansion is applicable to hadronic transitions of 
T(3S) into lower level bottomonia, we consider the possibility that T(35) has a 
D— wave component. This assumption leads to a natural explanation of the tttt spec- 
trum in T(3S) — > T(1S) nn. Consequences of this assumption on other hadronic and 
radiative transitions of T(3S) are also discussed in the same context. 
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1. Introduction 



It has been suggested recently by two of us ( S.C. and P.K. ) that the nn spectrum 
in T(35) — ► T(1S)tttt can be explained by including a D— wave amplitude for the 
dipion system JIJ. The most general amplitude for a spin-1 particle decaying into 
another spin-1 particle with the emission of two pions is given by 



M = A n e^e u 



(q 2 + B E,E 2 + C ml) + D { Vil p v + p^) , (1) 



in the lowest order in pion momenta expansion. Here, e and e are the polarization 
vectors of the initial and the final T's, p and p are the 4-momenta of two pions, 
q 2 = (p + p) 2 = m 2 ^ = s w , and Ex and E 2 are the energies of each pion in the 
rest frame of the initial T. Two sets of parameters give the best fit to the 
distribution in T(35) — > T(15 , )7T7r IJ. Various angular distributions of the decay 
products in e + e _ — > T(3S) — > T(1S)ti7t are predicted in Ref. |J, and these need to 
be verified by future experiments. 

However, the reason for D ^ in ([J) was not clearly discussed in Ref. Two 
possibilities were briefly mentioned : either a D— wave admixture in T(35) or, a 
breakdown of QCD multipole expansion for hadronic transitions of T(3S f ). It is our 
purpose to explore the first possibility in detail. Since QCD multipole expansion 
enables us to understand hadronic transitions between heavy quarkonia other than 
T(3S'), it is desirable to try to understand the amplitude (|l|) in the same framework. 
If this is possible, then other hadronic transitions of T(3S") can be studied in the same 



context. We note that theoretical predictions on hadronic transitions of T(3S) in the 
literature are not reliable, since they do not correctly describe T(3S) — > T(1S)ti7y. If 
our predictions are in serious contradiction with the experiments, then we may have 
to conclude that the QCD multipole expansion breaks down in case of T(3S). 

This work is organized as follows. In Section 2, the amplitude ([!]) is interpreted in 
the framework of QCD multipole expansion. It is found that results in Ref. ]l| can be 
readily obtained, once T(3S) is assumed to be a mixture of S— and D— waves with 
a mixing angle 4> '■ 

\T(3S)} = cos0 \3S) + sm(f)\D}. (2) 

Consequences of this assumption on other decays of T(3S) are then explored in detail. 
First of all, it turns out that the current upper limit on B(T(3S) — > T(1S) + rf) 
selects P2 from the two sets of parameters of Ref. |J. In Section 3, various radiative 
transitions of T(3S) are considered. There, a tight constraint on the D— wave mixing 
arises from electric dipole radiative transitions T(3S) — > Xbj(2P) +7. In the presence 
of a D— wave component in T(3S), some new and interesting radiative decays appear. 
It can affect the decay rate of T(3S*) — > rjb + 7, and allows the following cascade 
transitions : 

T(35) -1+ l D 2 -1+ h b (lP) -1+ Vb . (3) 

Besides these decays, T(3S') — > hb(lP) + 7r° and T(3S) — > hb(lP) nn are also inter- 
esting, and the D— wave contributions to these processes are considered in Section 4. 
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For these decays, we adopt the approach proposed by Voloshin 0, which correctly 
predicts the ratio of the charmonium 1 P 1 state decaying into J/ij) + 7r° and J /if) + 7T7T 
||. All of these decays reach a spin-singlet P— wave state, /if, (IP), that is hard to 
produce in the e + e~ annihilation. /if,(lP) can be a source of a spin-singlet S 1 — wave 
state (?/&) through electric dipole radiative transition, /if, (IP) — > r^f, + 7. Finally, our 
results are summarized in Section 5. 

In the following, the absolute decay rate or its lower/upper bound is derived 
for each decay process. It depends on the mixing angle <fi and quarkonium matrix 
elements of operators, r and r 2 , where r is the spatial separation of b and b. The 
matrix element of r 2 between T(35) and T(15) can be directly extracted from the 
spectrum and the absolute decay rate of T(35) — ► T(15) nir, and gives information 
that is independent of specific potential models. Absolute decay rates of T(35) — > 
Xt,j(2P) + 7 give useful information on the matrix element of r between T(3S) and 
Xf>(2P) and the mixing angle <fi. Other unknown quarkonium matrix elements will be 
fixed by the results from potential model calculations. We use m& = 4.8 GeV in this 
work. This induces some uncertainty less than ~ 20% in the numerical estimates of 
1/m 2 . Finally, some of our results in Sections 3 and 4 show explicit dependence on G 8 , 
the Green's function of the color octet bb states (defined in (0).) These results should 
be regarded as being order-of-magnitude estimates because of the approximation we 
will make about G$. 
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2. Hadronic transitions of T(35) into T(IS') 

Let us begin with the nir spectrum in T(3S) — > T(1S , )7T7t. In QCD multipole 
expansion, this process occurs through El — El multipole interaction, where the El 
interaction Hamiltonian of quarkonium with a gluon is [|J] 

H mt (El) = - l - gC n E?(0). (4) 

Here, g = (47ra s ) 1//2 is the SU(3) C gauge coupling constant, r is the relative position 
of the quark (6) and the antiquark (b), and £ a = t^ — tf is the difference of the SU(3) C 
color generators that act on the quark and the antiquark, respectively. When acting 
bewteen colorless states, £ a satisfies 

, 2 h 

( colorsinglet | ^ a ^ | colorsinglet ) = - 5 . (5) 

From @, the amplitude for i — > fun is given by 

M {i -> /tttt) = I ( / I n G 8 Tj I % ) ( tttt | 7ra s E?Ej \ ), (6) 

where G s is the Green's function for the color octet QQ states : 

Here k runs over color octet QQ states only. E and are the energies of the initial 
and the intermediate states. G$ is unknown due to our ignorance of quark confinement 
in QCD, and will be treated as a constant. Then, a lower bound on Gg can be derived 

IGsl 2 > 18 GeV~ 2 , 



using a sum rule on (nS \ r 2 \ IS) and the decay rate of T(2S) — > T(1S) irir 0. 
G 8 can be determined if the absolute decay rate of h c (lP) — > J/ip + ti° is known 
experimentally ||. 

Now we show that the D— term in naturally arises from the | AL| = 2 transition, 
3 Di — > 3 Si mi. We assume that the D— wave mixes with the initial quarkonium 
T(3S), since the final quarkonium T(IS') is the lowest level bottomonium and it is 
hard to imagine that it would contain any contamination of a D— wave. According to 
potential model calculations @, there are two D— wave levels below BB threshold, 
with m(lD) = 10.16 GeV and m(2D) = 10.44 GeV, respectively. Since T(2D) 
is closer to T(35) than T(1D), one might guess the D— wave in (|2]) to be T(2D). 
However, as discussed below (p^), a larger value of (IS* | r 2 | D) is desirable for 
mixing. Thus, T{1D) may be preferred because it has no node in the radial wave 
function. In this work, we do not address questions regarding the origin of the S— and 
D— wave mixing and which of the two D— wave levels enters in (|^). The discussion in 
this section is independent of such issues. In the next section on radiative transitions 
of T(3S), we consider both ID and 2D mixing. 

The angular part of the matrix elements between quarkonia can be easily per- 
formed, and we get 

(( 3 S 1 )' j | r k G 8 n | ( 3 5' 1 ) i ) = l - I SiS < 5 i3 S u , (8) 

\f^2i f 2 \ 

({ 3 S 1 ) j | r k G s r { \ { s D 1 ) i ) = — I S , D (S ik Sji + S u S jk - - 6 {j 5 kl J , (9) 
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where 

hj = (f I r k G 8 r k \i) = R f r k G 8 r k R { r 2 dr. (10) 

J 

The gluonic matrix element, {nit \ 7ia s E^E® | 0), can be calculated by considering 
(nn | a s G^p G va | 0) and the QCD scale anomaly. Detailed procedures are discussed 
in Refs. ||, [0] and The result can be summarized as 

(tt+tt- I a s Gp P G va | 0) (11) 
= A (q 2 + ml) + B (sw - # M<T T pi/ + g pcj Tp U - g pv T^), 

where 

V = PvPv + VvVp = 2 (Qui* - V")j ( 12 ) 

A = i^ = -i* (is) 

3/539' V ; 

5 = 5 «*Pg = i A, (14) 

and q — p+p , r — p — p . 

The A— term receives contribution from the QCD scale anomaly Q, while the 
B— term arises from the gluonic contribution to the energy momentum tensor of QCD 
||. The parameter A can be determined from the tttt spectrum in T(25) — > T(15) tttt 
H : A = 1.6 ~ 1.9. This is consistent with what we obtain below from the tttt 
spectrum in T(3S f ) — > T(IS') 7r7r. 

Using the information given above, one can calculate the S— and D— wave con- 
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tributions to T(35) — > T(lS)im : 



M(3S — > IS tttt) 



4tt 2 



x 



hs,is e • e 

51 

? 2 + ^-^ {(g°) 2 -(r°) 2 + g 2 + 2m 2 } 



M(kD -> IS tttt) 



a/2 A 
30 



(15) 



(16) 



M(T(3S) ->T(lS r ) tttt) =7W(35^ 15 tttt) cos0 + .M^D 15 tttt) sin0. (17) 
Note that the structure of the D— term in ([!]) comes from the first curly bracket in 



(pig), as mentioned at the beginning of this section. 



We fit the 7T7T spectrum in T(3S) — > T(15) nn using the above amplitude with 
three free parameters, IkDis sin0, /35,1s cos0 and A. The best fit is given by two 



sets of solutions (see Fig. 1) : 



lkD,lS 

hs,is 



tan0 = ±(2.4 ±0.5), 



A = (2.0 ±0.1) 



with x 2 /d-°-f- — H-2/7 (equivalent to 13.2 % C.L.). These correspond to two best 
fits (called PI and P2) obtained in Ref. |TJ using amplitude (W). More specifically, 
one can express (|T7|) in the form of ([!]) using ( fl5|) and (|16D , and find the value of D. 



It turns out that the upper and the lower signs in flilf) correspond to the parameter 
set P2 and PI (See table 1.) in Ref. [fjj], respectively. They can be distinguished by 
measuring various angular distributions of the final decay products as suggested in 



Ref. |J. Also, as discussed below in detail, the decay rate for T(3S) — > T(IS') ±77 can 



resolve this twofold ambiguity, and the parameter set P2 is preferred. The value of 
A = (2.0 ±0.1) obtained here is consistent with the A extracted from the 7T7t spectrum 
in T(25) -> T(15) 7T7T §. 

From the absolute decay rate of T(35) — > T(15) nn, we obtain the absolute values 
of I 1D) is sin and 135,15 cos : 

|/35,i5 cos 0| = \ (1S \ r G 8 r \ 3S)\ cos0 « 0.78 GeV~ 3 (19) 
|4 Djl 5 sin 0| = I (IS' I r G 8 r \ kD)\ sin0 « 1.92 GeV~ 3 (20) 



tan0 = 2.46 



1(15 


rG 8 r 


3S)\ 




r G 8 r 


1 ^)| 



(21) 



Up to now, we didn't care about G 8 . If we assume that G 8 is a constant, (|2TD becomes 



tan0 = 2.46 



1(15 


1 


35)1 


\(is\ 




1 kD)\ 



(22) 



Using the values of quarkonia matrix elements quoted in Ref. [0], 

|(15 I r 2 I 35)| = 0.3 GeV" 2 , 
|(15 I r 2 I kD)\ = 1.65 GeV" 2 , (23) 

we get ~ ± 24°. However, as discussed in Ref. 0, the matrix element | (15 | r 2 | 35) | 
can be much smaller than the above number, since the 35 state has two nodes which 
may lead to almost complete cancellation. Furthermore, the accuracy of the wave 
functions determined in potential models is about 10%. Therefore, the actual mixing 
angle may be much smaller than 24°. Also, because of our approximation on G 8 , 
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this kind of determination of is less reliable than that obtained in the next section 
from radiative decays, T(3S f ) — > Xbji^P)+l- In fact, too large a mixing angle (0) may 
result in severe discrepancy between the theoretical predictions and the experiments 
for these radiative decays. 

We can also consider the similar decay T(3S f ) — > T(2S) tttt. The best fit for = 0° 
(no D— wave mixing) yields ^/d.o.f. = 1.7/7. If the D— wave mixing is allowed, we 
get the best fit with y 2 /d.o.f. = 1.4/6. Therefore, it is difficult to tell which of the two 
is a better fit, and our assumption on the D— wave mixing cannot be tested clearly in 
this decay mode. Improved measurements of the 7iir spectrum in T(3S) — > T(2S) tttt 
are welcome. 

The message of this work can be put in the following way : a small mixture of 
a D—wave component in T(3S) can explain the tttt spectrum in T(35) — > T(IS') titt 
in the framework of QCD multipole expansion, if \ (1S \ r Gs r \ 3S)\ is much more 
suppressed compared to \ (1S \ r Gs r \ kD)\. 

One can also consider another transition, T(3S) — > T(IS') + f] (or n ), which oc- 
curs through interference of El and M2 interaction. The M2 interaction Hamiltonian 
is given by 

H int (M2) = ~ g Sj C nDiHm, (24) 
4m J 

— * 

where m is the mass of the quark, S is the total spin of the quark and the antiquark, 
and Di is the spatial component of the covariant derivative. 
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From (U) and fl24|), we can derive the amplitude for this transition and calculate 
the gluonic matrix element using the ?7a(1) anomaly in QCD M : 



M(?(3S) ^T(1S)+T]) 
= i (d k (r) | 7ca s E^H^ | 0)) m Q l e ij{ e, e\ I 3S ,is cos<f) - 4d,is sin. 
= y (J) /tt mg 1 e iife q ^- (p) k I 3S> is coscp- I k D,is sin. 
where Uj is defined in (|10|). Using the results (|19|) and (p0|) , we predict 



(25) 



r(T(35) -> T(LS) +77) 



58 eV (for P2), 
870 eV (for PI), 



or 0.2% or 3.6% in the branching ratio. Current upper limit on this decay mode is 
0.22% 0], which prefers the first set (P2) : 



B(T(3S) -> T{1S) + 77) = 0.2% (for P2), (26) 

which is close to the current upper limit. Therefore, twofold ambiguity encountered 
in Ref. [0 is lifted in the present work, and the parameter set P2 is preferred. Obser- 
vation of T(3S) — > T(1S) + 77 at the anticipated branching ratio would constitute one 
of the cleanest tests of our assumption : applicability of QCD multipole expansion to 
T(3S), and a small admixture of D— wave component in T(3S). 



3. Radiative transitions of T(35) 

In order to further check the D— wave mixing in T(3S), we consider electric dipole 
radiative transitions, T(3S) — > Xbji^P) +7 and T{3S) — > Xbj(^P) +7- In this section 
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and the following one, we assume the D— wave in T(3S) can be either ID or 2D state 
and consider both possibilities on the same footing. The transition rate of these 
decays is given by 



r(T(35) - Xb j(nP)+l) = — aQ 2 b uj 3 (2 J, + 1) \{nPj \ r | T(35))| 2 , (27) 



where 



(nP | r | T(35)) = (nP | r | 35) cos0 + V2 (nP \ r | fc£>) si 



smi 

(nP 1 | r | T(35)) = (nP | r | 35) cos0 - (nP | r | fcD) sin0, (2* 

V2 

(nP | J" | T(35)) = (nP \ r \ 3S) cos0 + — j= (nP \ r \ kD) sin0. 



The measured rates of Y(35) — > x&j(2P) + 7 are available JTT[] for J = 0, 1, 2. Anal- 
ysis of these decay rates yield a solution for the two variables (2P\r\3S) cos (ft and 
(2P\r\kD) sin0 : 

(2P | r | 35) cos0 = +(2.66 ± 0.16) GeV -1 , 

(2P | r | fc£>) sin0 = -(0.14 ±0.18) GeV" 1 . (29) 

To determine <fi from (p9f), we use the potential model calculations of ( / | r | i ) 
given in || : 

(IP | r | 35) = -0.023 GeV" 1 , 

(IP | r | IP) = -2.0 GeV" 1 , 

(IP | r | 2D) = -0.26 GeV" 1 
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(2P | r | 35) = +2.7 GeV~\ 



(30) 



(2P | r | ID) = +1.9 GeV~\ 
(2P | r | 2D) = -2.7 GeV -1 . 

Multiplicative relativistic correction factors to the matrix elements involving the 
S— wave have been calculated in Ref. |jl2| . The corrections depend on the state of 



Xbj(nP) into which the T(3S) decays : 



(2Pj | r | 35) : 



(IPj | r | 3S) : 



J = 2 
J= 1 
J = 
J = 2 
J= 1 
J = 



1.02, 

1.00, 

1.95. 

2.3, 

1.2, 

1.9. 



(31) 



From the second equation of (|29D and (|30|) , we get a value for the D— wave mixing 



angle as deduced from T(3S') decaying into Xbj(2P) : 



ho = -4° ± 6 C 
h w = +3° ± 4 C 



(32) 



The error in the angles is an estimate only. Moreover, the angles are seen to be 
consistent with zero. may also be estimated from the first equation of ( p9|) and 
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Eqs. ( pO|) and fl3ip . —12° < < +12° obtained this way is not very useful because 
the allowed range is large. 

A determination of the mixing angle is also possible by making use of experimental 
bounds on various combinations of branching ratios for T(3S) decaying into Xftj(l-P)- 



Form Ref. [|IT| we have the experimentally measured 



B(xto(lP)^T(lS)y) = 0.22 ±0.04 

B(xbi(lP) -> T(15)7) = 0.35 ±0.08 (33) 
B( Xb0 (lP) - T(15-)7) < 0.06. 

S(T(35') — > xw(l-P)7) is sensitive to the D— wave mixing angle and can be calculated 
from Eq. (P?D, ( p0|) and (|3T| ) knowing the total decay width. These branching ratios 



lead to a bound on <fi because one has to satisfy the following experimental relation 
P"B| (individual branching ratios are not available at this time) : 

F{<j>) = £ B(T(3S) - Xbj{lP)l) ■ B( XbJ (lP) - T(15) 7 ) 

J=l,2 

= (1.2 +^±0.09) x 10" 3 . (34) 

In Figs. 2 (a) and (b) we show F{4>) for mixing with \1D) and \2D) states respec- 
tively. The allowed region of <ft is larger in the case of 2D mixing if one demands 
the consistency between the Xwip-P) decays. We like to emphasize that the present 
experimental data is consistent with the assumption of a D— wave mixing. Based on 
our analysis, mixing with the 2D state is seen to be more plausible. 

13 



An estimate of the mixing angle (^) may now be used to predict the branching 
ratio of Xbo decaying into T(1S) using [13| 



£ B(T(3S) -> XW (1P) 7 ) • B(xu(lP) - T(15)7) 

J=0,l,2 

= (1.7 ±0.4 ±0.6) x 10~ 3 . (35) 



We finally write down expected branching ratios of T(3S f ) — > x&j(1D) + 7 f° r 
J = 0, 1, 2 assuming |1D) and |2D) mixing using the central values of Eq. ([32]): 

B(T(3S)^ X m(lP)+l) =1-1% (ID), 1.8% (2D), 

B(T(3S)^ Xb i(lP)+l) =2.0% (ID), 0.2% (2D), (36) 
D(T(35)^X6 2 (1D)±7) =0.3% (ID), 1.0% (2D), 

It is clear that a better determination of branching ratios of these radiative decays, 
T(3S) — > Xbj{nP) + 7 with n = 1, 2 can resolve 2D mixing from ID mixing, or vice 
versa. 

The spin-flip radiative transition T(35) — ► rji, + 7 is also affected by the D— wave 
component in T(3S f ). The decay rate is given by 



2 uj 3 

n 



r(T(35) ^ ^ ± 7) = - « Q — |ff , (37) 

3 mr 



where 

D = V2 (l^l jo(ujr/2) 1 35) cos0± j 2 {uJr/2) \kD) sin0, (38) 
where j n (x) is the n— th spherical Bessel function. In the long wavelength limit 
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(ur — > 0), we can approximate: 



j (x) = 1 - ^r, 
o 



-ttt hs,is COS + — J fcjD)ls sin. 
24 oU 



so that 

1 T a/9 1 

(39) 

assuming G 8 is a constant parameter, as in Ref. 0. Therefore, one can again use 
(|Hj|) and ( P0| ) to evaluate F, and calculate the decay rate from Q3"7| ) : 



4 (l8/G|(GeV~ 2 )J eV (for P2) 
75 (l8/G2(GeV- 2 )) eV (for PI), 



r(T(3S) - ^ + 7 ) = { _ V, „ , ,- 2 A , , ,, ™" ( 4 °) 



Note that (i) there is no 0— dependence left over, once we use the results (|T9|) and 
( p0|) , and (ii) this result is independent of which of the .D-wave actually mixes. If the 
parameter set PI were the correct one and G% not too large, this could open up a 
new option for the discovery of r] b in the e + e~ annihilation bypassing the intermediate 
stage involving fo&(lP). Unfortunately, the upper limit on B(T(3S) — > T(15) + rj) 
prefers the parameter set P2, for which B(T(3S) —>■% + '-/) < 1.6 x 10~ 4 . Therefore, 
this channel may compete with other possibilities discussed below in the search of T]b, 
only if G\ is not too large. 

The D— wave component in T(3S) can generate other interesting radiative tran- 
sitions : 

T(35) -1+ l 1 D 2 -1+ hb (lP) Vb . 
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The first chain is energetically allowed only for ID mixing. These decay rates can be 
readily obtained from the results of Ref. f[4" |. Omitting all details, the final results 



are given below : 

r(T(3S)^ 2+7 ) = f^^sinV (41) 

„ 3. 8 (^)y (**V eV, 

V 200 J v o- 1 / 

r(l 1 D 2 -> h 6 (lP) + 7) = ^ « Ql^ I (li 3 k I ID) | 2 « 30 keV, (42) 

15 

r(/i fe (lP) -> 776 + 7) « 40 keV. (43) 

This corresponds to the production of ~ 80 /i&(1P)'s in decays of 10 6 T(3S)'s for 
sin0 = ±0.1, or (p = ±6°, and if B{l l D 2 — ► h b {\P) + 7) = 50%. These are sensitive 
to the mixing angle 0, and may be useless for producing h b {lP) and r] b , if sin0 < 0.1. 

4. Hadronic transitions of T(35) into hb(lP) 

Finally, let us consider T(35) — > h b (lP) + X with X = 7r° or mi. This may 
serve as a source of the spin-singlet P— wave state, h b (lP), if its branching ratio is 
appreciable. In QCD multipole expansion, the above transitions are generated by the 
interference between PI and Ml interactions, where the Ml interaction Hamiltonian 

H int {Ml) = ^-gC A.tff (0), (44) 
2m 

— * 

with A being the difference of the spin operators for the quark and the antiquark. 
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The amplitude for T(35) — > hb(lP) + X is 



M(T(3S) - h b (lP)+X) = - ^- £i e'j ( X \ n a s E? H a k | 0) 



x( {h b (lP))j | r, G 8 A k + A k G 8 n | (T(35))i }, (45) 

where e and e are the polarization vectors of T(35) and hb(lP), respectively. 

The angular part of the quarkonium matrix elements can be performed as before, 
and we get 

( CPi)j I n G s A k | CS^ ) = i= (P P | r G 8 | R s ) 5 lk 5 jh (46) 

x/3 

( (^Oj I r, G 8 A fc | ( 3 L>0i ) = ^= (P P | r G 8 | R D ) 



x f<% 4z + 5u 6j k - - Sik Sji) • (47) 



Now, consider the case X = ir°, for which the matrix element of the gluonic 
operators are determined by Ua(1) anomaly and the mass difference between u and 
d quarks § || : 

(vr | vr a s E? H a k | 0) = 5 lk -f= (^LZp£\ U m *. = A 5 lk . (48) 

3V2 \m u + m d J 

Using the pion decay constant f n = 132 MeV and (m u — md)/(m u + m^) = 0.3, we 
get A = 1.7 x 1(T 3 GeV 3 . The amplitude for T(35) -> /i 6 (1P) + tt° becomes 

.M(T(3£) -> /i 6 (1P) + 7t°) = A J w e • e , (49) 

where 

4 

X ° 7S (IP | r | 3S) cos0 + V2 (IP | r | IP) sin0 , (50) 



9 nib 
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and the decay rate is 



r(T(35) - h b (lP) + n°) = ±- (A Ijt) 2 \pn\. (51) 

Z7T 



The expression in the bracket of fl5"0|) is the same as the first one in (pip. Therefore, 
T(35) — > xm(I-P) + 7 an d T(35) — > h h {lP) + 7r° are related with each other. The 
ratio of the decay rates of these two seemingly different decays is independent of 
quarkonium matrix elements of r or the mixing angle 0, and determines Gs- From 
(|7|) and (HD, we find 



T(T(3S) -> fefc(lP) +tt°) 27 / 4yg ggA 

r(T(35) -*xw(lP)+7) ~ 8naQ 2 b \1T 



2 

-0 \ 



3 (Gl (GeV- 2 - 



« 8.7 x 1(T 3 (^^^ J -j . (52) 

The absolute decay rate of T(3S) — * hb(lP) + 7T° is 

(For numerical estimates, we have used the values of (lP\r\3S) and (lP|r|lZ?) quoted 
in fl3U| ) with relativistic correction factors fl3"T| ) and central values of (52).) This 



amounts to the branching ratio greater than 4 x 1CT 5 for ID mixing, and 1.2 x 10~ 5 
for 2D mixing. Therefore, this decay may be the best for reaching hb(lP), and 
subsequently i] b through h b (lP) — > r\ h + 7. 

A similar decay, T(35) — > hb(lP)im, does not receive any contribution from the 
trace of the energy-momentum tensor in QCD, and is not enhanced over T(3S) — > 
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h b (IP) + vr°. From the general expression (|TT])-(0), we get 

( 7T+7T- | ix a s E? Hf I 0) = ~ X e ijk (E lP2k + E 2 p lk ), (54) 

where p^ = (Ei,px),p' IJ ' = (E 2 ,p 2 ) are the four-momenta of the pions. Then, the 
decay rate for T(3S) — > h b (lP)n + n^ is 

r(T(35) - ^(lP)7r+7r-) = ^ ^ |/ 27r | 2 , (55) 

where A = m(T(3S)) — m(h b (lP)), ip = 0.22 is the suppression factor of the phase 
space integral due to the pion mass 0, and 



4^3 G s 

J-2-k — 



(IP I r | 35) cos0 -= (IP | r | ID) sin. 

V2 



(56) 



9 m b 

Again, the expression in the square bracket is the same as the second equation of 



(HI), and this decay is related to T(3S) — > Xw(l-P) + 7 in the same way as in (|52|) . 
From (|BTD and (|55|) , we get 

r(T(3g) -> h b (lP) tttt) ^ f 0.2 (for ID) 

r(r(3S)^h b (lP)+n»)~\ 0.03 (for 2D) 1 J 

where we have used A = 2 and Eqs. (|30|), (|3~1~D and (0). This corresponds to ~ 10~ 5 
(for the ID mixing) and ~ 4 x 10~ 7 (for the 2D mixing) in the branching ratio of 
T(3S) — > h b (lP) ixix. The current upper limit is about two orders of magnitude above 
our prediction fllOfl . Our analysis shows that one has to look for T(3S) — > h b (lP) + 7r° 
(10~ 5 ~ 10 -4 in the branching ratio) rather than T(3S) — > h b (lP)7T7i to identify the 
spin-singlet P— wave bottomonium state, h b (lP), even if T(3S) has an admixture of 
a D— wave. 

19 



5. Conclusion 



In conclusion, we find that the irn spectrum in T(3S) — > T(1S) ttit can be ex- 
plained in a natural way in the framework of QCD multipole expansion by assum- 
ing the physical T(3S) state is an admixture of the S— and D— waves. The nil 
spectrum determines the mixing angle to be less than ~ 24°. (Decay rates for 
T(35) — ► xw(2-P) +7 give much tighter and more reliable value for <ft to be —(4° ±6°) 
for ID mixing, and +(3° ± 4°) for 2D mixing.) Effects of this D— wave component 
on other decays of T(3S) are discussed in detail. Twofold ambiguity encountered in 
Ref. I] is resolved by the upper limit on B(T(3S) — > T(1S) +rj), and the parameter 
set P2 is preferred for which 

B(T(3S) -> T(1S) + V ) = 0.2%. 

This is one of the cleanest tests of our assumptions on the D— wave mixing in T(3S) 
and applicability of QCD multipole expansion to hadronic transitions of T(3S). Con- 
sistency of our approach can be cross-checked by measuring the polar angle distirbu- 
tion ( with respect to the e + e~ beam direction ) of T(IS') in T(3S) — > T(15) nn as 
discussed in Ref. Further tests of our assumptions will be possible by checking 
our predictions (|36D by improving measurements of B(T(3S) — > Xbj{nP) +7). More 
informations on <fi and/or G 8 can be extracted by checking fllPD , (|4lD , (42|), (52) and 



(^3|). The most promising place to look for hb(lP) may be T(3S) — > hb(lP) +vr°, but 
the branching ratio is rather small, (|53|). (If sin0 « 0.1 and ID mixing is the correct 
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one, the cascades (3) may be useful.) Because of the small branching ratio of decay 
modes involving h b (lP) in the final state, it may be worth while to look at the direct 
transition T(3S f ) — > r/& + 7 for rjb search. 

In this work, we have derived several new results on the branching ratios of 
hadronic and radiative transitions of T(3S). If any of our predictions is in seri- 
ous contradiction with the experiments, then we may conclude that QCD multipole 
expansion breaks down for hadronic transitions of T(3S f ) to lower quarkonia because 
of the large radius of Y(3S) and its being close to the BB threshold. 
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Figure Caption 

Fig. 1 The best fit to the tttt spectrum in T(3S) — > T(IS') tttt using the amplitude 
(|17|). The results of this fit is given in (|i~8"D , for which x 2 /d.o.f. = 11.2/7. 

Fig. 2 Plots of F((p) of Eq. (|34]) assuming the mixing angle {(f)) to be a parameter. 
Fig. 2 (a) assumes a mixing with \1D), while Fig. 2 (b) assumes a mixing with \2D). 
The shaded region is allowed by Eqs. (|32"D and fl3~4"|). These may be used to obtain 
bounds on 0. 

Table 1. 

Two sets of parameters giving the best \ 2 fit to the tttt spectrum from Ref . |l| . 



Parameters 


Fit 1 (PI) 


Fit 2 (P2) 


A 


101.60 ±103.7 


366.34 ±61.94 


B 


-5.80 ±8.29 


-3.12 ±2.67 


C 


20.53 ±84.09 


4.33 ± 24.86 


D 


3.73 ±4.35 


-1.03 ±0.30 
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